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ABSTRACT 

Electron-positron pairs may be produced by either one or two photons in 
the presence of a strong magnetic field. In magnetized plasmas with 
temperatures kT ~ me 2 , both of these processes may be important and could be 
competitive. We calculate the rates of one-photon and two-photon pair 
production by photons with Maxwellian, thermal bremsstrahlung, thermal 
synchrotron and power law spectra, as a function of temperature or power law 
index and field strength. This allows a comparison of the two rates and a 
determination of the conditions under which each process will be a significant 
source of pairs in astrophysical plasmas. We find that for photon densities 
n < 10 25 cm“^ and magnetic field strengths B >_ 10^ G, one-photon pair 

'I * 

production dominates at kT > me 2 for a Maxwellian, at kT > 2 me 2 for a thermal 
bremsstrahlung spectrum, at all temperatures for a thermal synchrotron 
spectrum, and for power law spectra with indices s < 4. 
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I. INTRODUCTION 


Relativistic and trans-relativi Stic plasmas, where temperatures approach 

and exceed the electron rest mass energy, have been of increasing interest in 

astrophysics. The existence of such plasmas has been inferred from the study 

of sources such as gamma-ray bursts, where observed spectra require 

2 

temperatures of order me , and active galaxies, whose gamma-ray emission 

2 

implies temperatures in excess of me . These plasmas may have significant 
concentrations of electron-positron pairs, whose presence could strongly 
affect the spectrum, energetics and dynamics of the source regions. Study of 
the microphysics of pair production and annihilation is therefore of central 
importance to understanding the properties of hot plasmas. There has recently 
been a substantial amount of work on the processes governing the behavior of 
non-magnetized and weakly magnetized relativistic thermal plasmas. Rates of 
photon-photon pair production at relativistic temperatures have been 
calculated by Weaver (1976) and Stoeger (1977), and rates for particle pair 
production processes have been studied in the relativistic limit by Lightman 
(1982). Pair annihilation in hot thermal plasmas has been investigated by 
Ramaty and M£sz3ros (1981), Zdziarski (1981), and Svensson (1982a). Using 
results of these studies, Lightman (1982), Svensson (1982b), Stepney (1983), 
Takahara and Kusunose (1983) and McKinley and Ramaty (1983) have determined 
the steady state equilibria of a thermal plasma, where pair production would 
balance annihilation. 

The properties of trans-relativi Stic plasmas in strong magnetic fields 
are virtually unexplored, but have potentially interesting applications to 
gamma-ray burst source models. In comparison, the characteristics of strongly 
magnetized plasmas at lower temperatures (10 - 100 keV) have been extensively 
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studied in connection with models for X-ray pulsars. The coincidence of 
fields around 10 12 G with plasma temperatures around me 2 introduces first 
order processes, one-photon (photon-magnetic field) pair production and 
annihilation, which may be as important as the two-photon processes which 
operate in both magnetized and non - magnetized plasmas. In addition, 
synchrotron emission will be far more important as a source of pair-producing 
photons than bremsstrahlung or annihilation radiation. As a result, the 
behavior of strongly magnetized pair plasmas could be quite different from 
what has thus far been determined for non-magnetized pair plasmas. In view of 
the evidence for gamma-ray burst source regions being near the surfaces of 
strongly magnetized neutron stars (Katz 1982, Lamb 1981), it seems worthwhile 
to evaluate quantitatively the relative importance of the one-photon processes 
in a trans-relati vistic plasma. 

This paper calculates and compares the one-photon and two-photon pair 
production rates for various thermal and power law photon distributions in the 
presence of a strong (> 10 12 G) magnetic field. We consider pair production 
by photons with Maxwellian, thermal bremsstrahlung, thermal synchrotron, and 
power law energy spectra, with preliminary results having been given for the 
Maxwellian case in Burns and Harding (1983). With the exception of the two- 
photon rates for Maxwellian and bremsstrahlung photon distributions, none of 
these have previously been calculated. Guilbert (1983) has considered one- 
photon pair production by a power law photon spectrum, but did not give a 
calculation of the rate. In Section II, we present an outline of the 
calculations and a review of the less familiar one-photon pair production 
process. Sections III and IV give, respectively, the results for the one- 
photon and two-photon rates as well as the energy distributions of the created 
electron-positron pairs. We also to provide approximate analytic expressions 
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for the one-photon rates. In Section V, we look at the ratios of the two 
processes as a function of magnetic field strength and either temperature or 
power law index of the photon distribution. Section VI discusses the 
implications of the-results for self-consistent models of strongly magnetized 
pair plasmas and for models of gamma-ray bursts. 

II. OUTLINE OF CALCULATION 

The total rate of pair production by photon processes in a plasma is an 
integration of the cross section for the various processes over the photon 
distribution function. If the photons involved in the pair production are 
themselves produced by the particles of the plasma, the photon distribution in 
a self-consistent model would depend on the radiation mechanisms, the energy 
distribution of the particles, and the optical depth of the source region. 

The present calculation does not attempt self-consistency, but starts from 
assumed photon distributions which could plausibly be produced by mildly 
relativistic plasmas. This will give reasonable estimates for both the total 
rate of e+e- pair production and the relative importance of the one- and two- 
photon processes as a function of parameters which can be related to the 
properties of the plasma. 

In a strong magnetic field, the energies of the electron and positron 
perpendicular to B are quantized and the pair can therefore only be created in 
discrete states. When the photon center-of-mass energies are not large 
compared to the spacing of these states, the pair production rates are 
significantly influenced, primarily near threshold, by quantum effects. 

While an expression exists for the one-photon pair production rate near 
threshold (Daugherty and Harding 1983), none exists for the two-photon rate in 
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a strong magnetic field. For the sake of consistency in comparing the two 
rates, we will neglect the effect of the discreteness of the pair states, 
althouyh we will present additional calculations of the one-photon rates with 
quantum effects included. 

Throughout the rest of the paper, we will use the dimensionless 
variables: E for photon energy in units of the electron rest mass me 2 , E + for 
electron/positron energy in me 2 , 8' = B/B cr for the magnetic field strength 
in units of the critical field strength, B cr = m 2 c 3 /eh * 4.414 X 10 13 G, and 
temperature = kT/mc 2 . 

a) One-photon Pair Production 

The creation of e+e- pairs by single photons is possible in the presence 
of a magnetic field, because the field is able to participate in the momentum 
transfer, allowing momentum to be conserved (for a complete discussion of the 
kinematics, see Daugherty and Harding 1983). The probability of one-photon 
pair production, however, is negligibly small unless the quantity, x = 1/2 E 
B'sinQ > 0.1, where E is the photon energy, and 0 is the angle between the 
photon direction and B. The photon must also be able to supply at least the 
rest mass of the pair in the frame where 9 = tt/2, and must therefore have an 
energy in excess of the threshold, 2mc 2 /sin0. 

The polarization-averaged pair production attenuation coefficient for a 
photon of energy E propagating at an angle 0 to a constant, homogeneous field 
of strength 8' is (Toll 1952, Erber 1966) 

a Tg = j “ B' si n0 T (x) (1) 


T (x) - 4.74 x' 1/3 Ai 2 (x" 2/3 ) 


(la) 
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- 0.46 exp C-4/3x], X<<1 

(lb) 

- 0.60 x' 1/3 , X»1 

(lc) 


wnere Ai is the Airy function, * is the electron Compton wavelength and a is 
the fine structure constant. Eqn (1) is the asymptotic expression for the 
attenuation coefficient in the limit where the number of kinematical ly allowed 
states for the e+e- pair is large. Daugherty and Harding (1983) have found' 
that the change in the attenuation coefficient due to quantum effects in the 
case where the number of allowed pair states is not large can be approximated 
by making the substitution: 

X - x/F,' 

F = l + 0.42 (E sine)' 2,7 

in Eqn (1). The largest values of F, and thus the most significant changes in 
a Y(i , occur near threshold, where E sine = 2. The pair production rate "per 
photon" for a photon distribution, 

n Y (E, a) dEdfi = f(E, a) dEd a (2) 


wi 11 bs 




= ~ /d fl / n v (E. n) r y B (B* , E, Q) 


dE 


'mi n 


(3) 


where r Y g = c a T g is the rate for an individual photon and E m i n = 2/sine is 
threshold. The photon distribution is normalized to the total number density. 



8 


n Q , so that < is determined by 


If n y (E, a) dEdfl = n Q (4) 

The energy distribution of the created pairs can be found from the 
differential single photon rate (Daugherty and Harding 1983): 


dr YB (£ > 

cle 


si n8 


ctcB * /3~ [2+e (1-e)] 
< THttx E( 1 - e ) 


K 2/3 


[ lxcn7eT ] 


(5) 


where e = E + /E is the fractional energy of one member of the pair (note that 
the distribution is symmetric about e = 0.5 because the electron and positron 
are identical). The differential positron (electron) production rate for a 
distribution of photons is then. 


dR 


y B < E + ) 

3F 


oo 

= 2J- / ® / " v (E, n) 


dr YB (e > 

de 


de 


dE 


( 6 ) 


'mi n 


where E min = E + + 1/si nQ . 


b) Two-photon Pair Production 

Although two-photon pair production is a process which occurs in free 

space, it will be modified in the presence of a strong magnetic field to the 

extent that the quantization of pair states discussed above is important 

(Wunner, Herold and Ruder 1983). Daugherty and Bussard (1980) have derived 

the cross section for the inverse process, two-photon annihilation of pairs in 

a magnetic field, and find that it does not depart significantly from the 

1 3 

free-space cross section until B > 10 G. However, their calculation assumes 
that the pairs are in the ground state when they annihilate. Since this 



9 


initial state is a small part of the phase space which would be kinematically 
available to the pairs in the final state of pair production, we are not able 
from their results to estimate the general quantum behavior of the two-photon 
pair production rate. Since the cross section for photon-photon pair 
production in a strong magnetic field has not yet been derived we use the 
invariant free-space cross section (Berestetski i , Lifschitz and Pitaevskii 
1971): 


o (E. E\ <!>) - l°-l {(t 2 + 4t-8) In [ — — ] - (t+4) /t(t-4) } 

£ /t - /TT7T 

(7) 

where x (E, E ' , <i/) = 2EE' (l-cos40> o-j- is the Thomson cross section, and is 
the angle between the propagation vectors of the two photons. The pair 
production rate for a distribution of photons requires an integration of the 
cross section over both photon energies, E and E‘, and the solid angles for 
their propagation di rections: 


QO 00 


-p = | / da / dfl' / dE / dE' n v (E, n) n Y (E 1 , O') ( l-cos*) a ( E,E 1 , «) 


El, 


min n 


( 8 ) 


where E' min (E, <i>) = 2/[E(l-cos4>)] 


(9) 


is thresnold energy for the primed photon. In the case of an isotropic photon 
distribution, the above angular integrations can be reduced to an integration 
over just one angle, which is usually chosen to be 4', the angle between the 
photons (see Weaver 1976). 


c) Photon Distributions 

We calculate the rate of pair production for both processes using four 
different photon distributions in energy and angle: Maxwel 1 -Bol tzmann (or Wien 
spectrum) (MB), thermal bremsstrahlung (with gaunt factors ignored) (TB), 
thermal synchrotron ( TS ), and power law (P). In terms of the dimensionless 
variables T*, E and B', defined above, we take the following expressions for 


the function f(E,fl) in Eqn (2): 

f MB (E, a) = E 2 exp (-E/TJ (10a) 

f TB (E, 51) = j exp (-E/TJ (10b) 

f TS (E, fl) = exp [-( ^ — -— ) 1/3 ] (10c) 

T; B'sine 

f P (E, n) = E' S , s > 2 (10d) 


The expression for TS is the optically thin synchrotron spectrum from mildly 
relativistic electrons with a Maxwellian distribution at temperature T^ < 1 
derived by Petrosian (1981) in the limit E.>> T ^ B 1 . The dependence on 8, 
the angle between the emitted photon and the magnetic field, makes the TS 
distribution anisotropic, with a higher density of photons propagating 
perpendi cul ar to the field. The other three distributions are assumed to be 
isotropic, even though thermal bremsstrahlung in a magnetic field should be 
anisotropic to a certain degree (see Nagel and Ventura 1983). From Eqn (4), 
the normalization parameter < for the distributions in Eqn (10) are. 
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TB = Ej (£ l /T.) 

(lib) 

.TS n T 2 nl 
~ J '* 8 

(lie) 

c (1-S) 


1 

jJ l/> 

II 

Ql 

(lid) 


where E^ is an arbitrary lower energy cutoff in the the bremsstrahlung and 

power law spectra to make the integration finite, and E L is the exponential 

integral. The lower limit of E=0 has been used to normalize the TS spectrum 

even though emission is expected to fall off below the first harmonic, E=B 1 . 

Thermal broadening at T ~ 1 , however, will cause the effective lower energy 

★ 

cutoff to be E << B‘, and in any case, the cutoff is always below pair 
production thresholds. 

III. ONE-PHOTON RATE 

a) Analytic Expressions 

In the limit, x << 1 and E >> 2, it is possible to perforin the 
integration over the cross section to obtain approximate expressions for the 
one-photon pair production rate for the four photon distribution functions. 
Using the limiting expression (lb) for the magnetic single photon pair 
production rate, we can rewrite Eqn (3) as, 

i 

r yB * tP T- b ' /* d8 s'" 28 /’ f < E - 9 > t- WsTnff J (12) 

0 2 

sTnF 


Since all of the photon distributions are decreasing functions of energy and 
the pair production rate is sharply increasing with energy, the E integration 
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can be performed by the method of steepest descents. The integrand has saddle 
points at the energies: 


E 


E 


E 


E 


8T* 


MB _ 

o ' v 3B'sin0 


(lll'cinfl) 


1/2 


TB _ , 8T * x 1/2 
o " l 3B r sin0 j 

TS _ „ 8T * \ 1/ 2 

o ' ^TI r "s’Tril' 


p - 8 . > ? 
o ' fs B TTnF’ 5 * c 


(13a) 

(13b) 


(13c) 


(13d) 


for the different photon distributions. These are the energies at which most 
of the photons pair produce. The major contribution to the integral lies in a 
relatively small region about the saddle point and there is negligible 
contribution elsewhere. Since the use of the steepest descents method 
requires extension of the lower limit of the integral to 0, the result will be 
most accurate in cases where E 0 >> E ml - n (ie., when the saddle point lies well 
above threshold). The 9 integration can then be performed by noting that the 
integrand also peaks sharply about sin0 = 1, and by making an expansion about 
this point, retaining only the linear terms. The resulting pair production 
rates for the four photon distributions are: 


R^ 8 = 1.11 exp [-2 (3^-) 1/2 ] 


YB 


rTB . .313 £ CE ^)]- 1 B ’ 3/2 exp [-2( 8 ) 1/2 ] 


■3 = 1,65 ^r* 


8 i i/2 i 


(14a) 

(14b) 


exp [-2( 


( 14 c) 



13 


*S B « .511 


ac 


-S-l 


(^) S - 1 


r(|+D 




(s-l)e 


■ s B lS 


( 14d) 


It is an interesting coincidence that this calculation gives the same 
exponential dependence in the TS case as in the MB and TB cases. Zheleznyakov 
(1982) has performed a similar computation of the one-photon pair production 
rate for a thermal bremsstrahlung photon distribution, except that his 
expression was not integrated over angle. The approximations we have made 
here', x << 1 and E >> 2, with the expressions for E 0 from Eqn (13) substituted 
for E, give the following regions of validity for the above expressions: 


MB.TB: (T*B') 1/2 « 1, T* > | B' (15a) 

TS: (| T*B* ) 1/2 « 1, T*> | B' (15b) 


P: s > ^, SB' < ^ (15c) 

The above expressions for R can be modified to take quantum effects 
into account by considering the second term introduced in the exponential in 
Eqn (lb) by the substitution x x/F as small. This term then becomes a 
slowly varying multiplicative factor in the steepest descents analysis which 
is taken outside the integral and evaluated at E=E 0 and sine * 1. The result 
is a mulitplication of the expressions in Eqn (14) by the factor. 


MB, TB: 

exp 

[-1.186 B' * 85 /T 1 * 85 ] 
* 

(16a) 

TS: 

exp 

[-.09 B' * 85 /T 1 * 85 ] 

(16b) 


★ 
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P: exp [-.193 B' 2 * 7 s 3 * 7 ] (16c) 

These factors give a good approximation to quantum effects in the pair 
production rate within the regions of validity estimated by Eqn (15). 

However, the most significant changes in the rate will occur for combinations 

of B' and or s ou tside the region where the analytic expressions are 
accurate. 

b) Numerical Results 

According to Eqn (15), at low temperatures or in magnetic fields 
approaching the critical value, where most of the pai r-producing photons have 
eneryies near threshold, the analytic expressions and approximation to quantum 
effects derived above are not accurate. We have therefore also numerically 
integrated Eqn (3), explicitly including the threshold photon energy and the 
full expression [Eqns (l),(la)] for the asymptotic rate, which does not 
require x to be small. The values for the saddle points [Eqn (13)] and the 
second derivatives from the steepest descents analysis were used as an 
estimate of the peak and width of the integrand to determine the numerical 
limits of integration. Figures la- Id present the numerical results for the 
pair production rate "per photon", both with and without quantum effects 
included, as a function of T^ (or power law index s) and B' along with the 
dashed analytic curves from Eqn (14) for the same parameters. Since 
expressing the rate "per photon" requires normalizing to the total number of 
photons in the spectrum, the TB and P rates are dependent on the additional 
free parameter E^ which cuts off the spectrum at some low value of the 
energy. We will thus plot R^g = R^g k for these cases. However, if E[_ << T* 
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in the TB case or E|_ - 2 in the P case, the TB and P one-photon rates can be 
compared directly with the MB and TS rates which are exactly normalizable. It 
is most convenient to discuss the thermal and non-thermal distributions 
separately, since the parameterizations are different. 

The pair production rates for the thermal photon distributions all have a 
similar behavior, indicated by the same exponential dependence in the analytic 
expressions. This exponential dependence requires the product T^B' to be 
reasonably large (> .01) before pair production becomes important. It also 
means that for temperatures below a few times me 2 , small changes in the 
parameters can cause order-of-magnitude changes in the rate. The rates for 
the MB distribution tend to be larger than those for TB because there are more 
photons above 1 MeV at a given temperature in the case of MB. The TS rates 
are larger than either of the others because there are more photons in the TS 
distribution which propagate perpendicular to the field, where the probability 
for pair production is highest. The TS distribution therefore picks out the 
most favorable phase space for pair production. The analytic curves match the 
numerical ones reasonably well over the parameter space where they were 
estimated to be accurate [Eqn (15)], but they overestimate the rates at low 
temperatures and high field strengths, where energies are near threshold, and 
at very high temperatures, where the single-photon rate is no longer 
increasing exponentially with energy. The quantum effect of discrete pair 
states is strongest for low temperatures, where the pair producing photons are 
near threshold, and tends to decrease the rate of pair production (by limiting 
the phase space allowed in the final state). The effect is not as strong in 
the TS case, because the value of E 0 in Eqn (13c) indicates that the pair 
producing photons have energies higher above threshold. 

In the case of the P distribution, the temperature is replaced by the 
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power law index as the parameter determining the number of photons above pair 
production threshold. The rates for the P distribution are generally power 
laws in B' with the same index s (only negative) as the photon distribution. 
Quantum effects are important here only at the highest field strengths and for 
s > 3, where pair producing photons have lower energies. 

c) Pair Energy Distributions 

The differential one-photon rate which gives the energy distribution of 
the created pairs is too complicated to evaluate analytically, even in the 
limit x << 1. so we evaluate Eqn (6) for several cases by numerical 
integration. Figures 2 and 3 show the results for the MB and TS distributions 
at two temperatures, 0.1 me 2 and me 2 . Figure 4 shows the differential rates 
for the P distribution with two values of power law index. For nearly all 
values of T* and B', each member of the pair will receive half of the photon 
energy (ie., the differential rate in Eqn [5] for a single photon is sharply 
peaked at e = 0.5). Therefore, the pair energy distribution has its maximum 
at around half of E 0 , the most probable pai r-produci ng photon energy. The 
width of the distribution is the dispersion in photon energy convolved with . 
the dispersion of positron (electron) energy about e = 0.5. The TB 
distribution gives results which are very similar to MB, since they share the 
same values for E 0 and for the dispersion of pai r-producing photon energies. 
The TS spectrum produces pairs at higher energies for the same values of T^ 
and 8’, since the value of E 0 is larger than for the MB or TB cases. At low 
temperatures and higher field strengths, most of the pairs are produced "cold" 
(ie., with only their rest mass). Photons in a power law spectrum produce 
pairs with a much wider range of energies, the width being inversely 
proportional to power law index, but there is still a well defined maximum in 
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the distribution. 

It is interesting to note that in a number of cases, most of the pairs 
are produced with energies considerably greater than the average energy of the 
particles in the plasma. The average energy in a Maxwellian distribution of 
electrons either in equilibrium with photons in the MB distribution at kT or 
which are radiating photons in the TS spectrum is 3kT. The curves shown in 
Fiyure 2 indicate that pairs are being produced with energies as high as 20- 
30kT in the case of the low field strengths. The rate of producing these high 
eneryy pairs, however, is very low. 

IV. TWO-PHOTON RATES 


Deriving analytic expressions for the rate of photon-photon pair 
production by means of the steepest descents method is not accurate enough to 
give useful results. The approximation which must be made for the cross 
section depends on E 1 ^ 2 , a much slower dependence than the one-photon cross 
section which is an exponential. This produces an overlap with the 
distribution function which is neither sharply peaked nor sufficiently 
symmetric about the saddle point for the Gaussian approximation made by the 
steepest descents method. We are therefore limited to a numerical evaluation 
of the two-photon rate. 

In the case of the isotropic photon distributions, we need only perform 
integrals over the photon energies and the angle between their propagation 
directions. Results for the MB and TB spectra are plotted together in Figure 
5a, although could shift up or down relative to R^“ depending on the low 
energy cutoff in the spectrum which determines the normalization. However, in 
the limit E|_ << T*, the dependence should be weak. The curve 
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for reproduces Weaver's (1976) result, except that he computes the 
quantity <av> = 2 R /n Q , which differs by a factor of two. In order to 
directly compare the one-photon and two-photon rates we must compute the rate 
of pair production events per photon, which is 1/2 of the total rate of events 
in the case of the two-photon process. 

The rates for the TS spectrum are plotted in Figure 5b for three field 
strengths. Even though we have used the free-space cross section, the two- 
photon rates are now dependent on B* through the dependence on B'sine in the 
TS distribution function. It is. therefore necessary to integrate over three 
angles in Eqn (8), which we have chosen to be 9, 9‘, and $ = <t>‘ - <j>, in the 
coordinate frame where B lies along the z axis. The angle t is then related 
to 9, b', and $ through the law of cosines for spherical triangles: 
cos*|j = cos 0 cos0‘ + sine sinQ' cos<j>. The rate of two-photon pair production 
from a TS spectrum is significantly less than from the MB or TB spectra, 
because of the anisotropy of the distribution function. The departure of the 
TS rate from the MB rate is largest for low temperatures and low fields, where 
the beaming of photons in the TS spectrum is most pronounced. 

The rate for the P case is plotted in Figure 5c as a function of power 
law index. Like the TB case, the two-photon rate for the power law 
distribution is dependent on the low energy cutoff in the spectrum. Except 
for some dependence on s in <, the two-photon rate decreases with increasing 
power law index because the steeper spectra have fewer photons at energies 
around the maximum in the cross section. 



V. RATIO OF ONE-PHOTON TO TWO-PHOTON PAIR PRODUCTION RATES 


The ratio of rates for one-photon and two-photon pair production "per 
photon" can be computed for each value of field strength and either 
temperature or power law index. However, the total photon density, n Q , is 
introduced as an additional parameter, since the relative importance of the 
two processes will depend on the relative densities of real photons and 
virtual photons of the magnetic field. Figures 6a-6d show the ratio of the 
one-photon to two-photon rates as a function of kT and B' (or s and B 1 in the 
P case). On the left-hand vertical scale is the quantity n Q R^B^yy’ which 
shows the photon density at which the two rates are equal. On the right-hand 
side is the actual ratio of the rates for a fixed photon density of n Q = 10 25 
cm -3 . For the thermal distributions, we also plot the blackbody or thermal 
equilibrium photon density, n BB , which is an upper limit to the photon density 
at kT. Therefore, in the region of Figures 6a-6c above n BB , the one-photon 
rate always exceeds the two-photon rate. 

For the thermal photon distributions, temperatures below 1 MeV are of 
most interest to models of gamma-ray bursts. In this region, small changes in 
the magnetic field strength can produce order-of-magni tude changes in the 
ratio of the two rates. For example, for the MB distribution and a photon 
density of n Q = I0 25 cm -3 , the two-photon process dominates up to T* - .6 for 
B 12 = 2, whereas for B 12 = 4.4 the one-photon process dominates at all 
temperatures. The ratio of one- to two-photon rates is significantly larger 
at all field strengths and temperatures for the thermal synchrotron 
spectrum. In this case, one-photon pair production is the dominant mechanism 

12 oc o 

at all temperatures for B > 10 G, and n 0 < lO* 13 cm -0 . 

When the pairs are in thermal equilibrium with the photons at temperature 
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T ~ 8', most of the pairs will be in the lowest Landau state and we can 
★ 

compare our results for the ratio of one-to two-photon rates with those of 

Daugherty and Bussard (1980) for the inverse processes. They find that the 

one-photon and two-photon annihilation rates for pairs in the ground state are 

equal at 8 = 10* 3 G. Using Fig. 6a for a Maxwellian photon distribution, we 

find that at the thermal equlibrium photon density, n^, one-photon and two- 

photon pair production rates are also equal around B = IQ 13 G under the 

constraint that T = B' . 

★ 

VI. DISCUSSION 

It is of interest to examine the implications of the above calculations 
of pair production in mildly relativistic plasmas for models of gamma-ray 
burst source regions and spectra. The observed spectra of many bursts show 
line features at energies around 400 keV which could be red-shifted .511 MeV 
photons from positron annihilation. Fits to the continuum yield temperatures 
of 100 keV to 1 MeV for thermal bremsstrahlung spectra (Mazets et. al. 1981) 
and temperatures of 100 - 500 keV and B ~ 10 12 G for thermal synchrotron 
spectra (Liang 1982). These properties suggest that creation of electron- 
positron pairs must be taking place in the source regions of these objects. 
Two-photon pair production has been considered as the source of pairs in most 
models for bursts (Cavallo and Rees 1978; Ramaty, Lingenfelter and Bussard 
1981), but if strong magnetic fields are present our calculations indicate 
that one-photon pair production could be at least as important. The observed 
ratio of energy in the annihilation line to energy in the continuum ( about 
0.1) requires at least this fraction of energy of the plasma to be in pairs. 
The contribution of pairs from magnetic pair production of bremsstrahlung 
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photons in gamma-ray bursts was estimated by Zheleznyakov (1982), who found 
that this process yives the observed 1 ine-to-continuum ratio provided that the 
source region is optically thick to pair production, that all the pairs 
annihilate, and that continuum radiation from the pairs is not important. 
However, in neutron star magnetic fields, the synchrotron loss rates are much 
higher than bremsstrahlung loss rates, so that in estimating the rate of pair 
production, a thermal synchrotron spectrum should be used. Since the 
corresponding rate of pair production from synchrotron radiation is much 
higher at a given temperature, an optically thick source will give too large a 
1 i ne-to-conti nuum ratio. One is led to conclude that gamma-ray burst source 
regions, if they have strong magnetic fields, are not optically thick to pair 
production . 

Clearly, a more sophisticated approach is needed in order to model gamma- 
ray burst spectra. To self-consistently determine the concentration of pairs 
in the plasma and the rate of pair annihilation, it is necessary to include 
all important radiation and pair production processes, and optical depth 
effects, both from pair production attenuation and scattering. In a thermal 
pair plasma, ie. one where the pairs have a Maxwellian distribution, the pair 
density may be found from a steady state solution which balances pair 
production and annihilation. A simplified model of steady state pair 
equilibrium in a strongly magnetized plasma, where the only processes included 
are one-photon pair production by synchrotron photons and two-photon 
annihilation, is considered by Harding (1984). This calculation gives 
relatively high pair densities, ie, near or above thermal equilibrium 
densities, for values of field strength and temperature estimated for gamma- 
ray bursts, unless the source size is much smaller than a neutron star polar 
cap radius. Not included, however, were Compton scattering of photons by the 
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pairs and optical depth of the photons from pair production. 

The optical depth of the sources to pair production is also an 

observational consideration. Photons above threshold would be severely 

attenuated in sources having high pair production optical depths. The 

observed spectra do not show strong attenuation above 1 MeV; in fact, a 

number of burst spectra have significant flux up to ~ 10 MeV and are best 

fitted by thermal synchrotron or power law models (Nolan et. al. 1984). These 

photons in the high energy part of the spectrum must be produced in regions of 

low optical depth. If thermal synchrotron models are to explain the continuum 

pnoton production then parameters derived from fits to observed spectra must 

give a low optical depth for pair production. The results presented in this 

paper can be used to obtain an average optical depth for photons in the 

spectrum. The optical depth for one-photon pair production is <t D > ~ R _ i/c 

where i is the path length of the photon through the strong field region. 

From Figures lc and Id, <t > - ( .3 - 10 9 ) l , where a is in units of 10 6 cm. 

Yd do 

a typical size for the strong field region near a neutron star. It appears 
that the observed photons with energies above 1 MeV cannot be produced within 
fields >> 10 12 G. Even if the source of observed photons is in a low field 
region, optical depths for two-photon pair production could attenuate photons 
above ~ .511 MeV. The two-photon optical depth <x^> a (R /n Q ) Rn 0 / C » 
depend on both source size R and photon density n 0 . Using a typical burst 
luminosity, L s 10 37 (d/100 pc) 2 erg s" 1 where d is the distance to the 
source, we can estimate the photon density in the source as, 

= > = 10 23 cm 3 (d/100 pc) 2 

o ~ <£>Ac " T* Aio 

where Aio is the area of the emitting region in units of 10 10 cm 2 . This gives 
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for the pair production optical depth at T* = 1 of a bremsstrahlung spectrum 
(using Fig. 5a) : 

<t > = 10 4 -r — (d/100 pc) 2 . 

IT Ml o 

So, unless R/Aio < 10 4 cm. or the sources are very nearby, two-phot on pair 
production will also attenuate high energy photons. This argument has in fact 
been used to place distance limits on gamma-ray burst sources (Schmidt 
1978). However Figure 5b suggests one possible solution to the problem. At 
temperatures below - .3 me 2 , the two-photon rate for a thermal synchrotron 
spectrum decreases with decreasing field strength significantly below the 
values for an isotropic distribution. The resulting pair production optical 
depths are 

<T yy >TS = (10 ' 9 * 10 ’ 5) Al^ (d/10 ° pc)2> 

If the photons in burst sources are produced by thermal synchrotron radiation 
and the fields are not high enough to cause one-photon pair production, the 
photons in the > 1 MeV range may be able to escape without strict requirements 
on either source geometry or distance. 

We would like to thank J. K. Daugherty, F. C. Jones, D, Kazanas, J. M. 
McKinley and R. Ramaty for valuable discussions and comments. We also thank 
C. Jones for programming assistance. 
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FIGURE CAPTIONS 


Figure la - One-photon (magnetic) pair production rate per photon at 
temperature kT/mc 2 for a Maxwellian photon distribution. Curves are labeled 
with different magnetic field strengths in units of 10 12 Gauss. Solid lines 
are numerical calculations, dashed lines are the analytic results, both 
without quantum effects; dot-dashed lines are numerical calculations including 
quantum effects. 

Figure lb - Same as Fig. la but for a thermal bremsstrahlung spectrum. 

R D = R u *. 
yB yB 

Figure lc - Same as Fig. la but for a thermal synchrotron spectrum. 

Figure Id - One-photon pair production rate per photon versus magnetic field 
in units of the critical field for a power law photon spectrum. Curves are 
labeled with power law index, s. 

Figure 2 - Spectrum of pairs produced by the one-photon process from a 
Maxwellian distribution for two different temperatures and magnetic field 
strengths. The vertical scale is the differential rate of pair production 
normalized to the total rate, and the horizontal scale is pair energy in 
me 2 . Curves are labeled with values of magnetic field strength in units of 
10 12 Gauss. 

Figure 3 - Same as Fig. 2 but for a thermal synchrotron spectrum. 
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Figure 4 - Same as Fig. 2 but for a power law spectrum, for two values of 
power law index, s. 

Figure 5a - Two-photon pair production rate versus kT/mc 2 for 
Maxwellian, and thermal bremsstrahlung, photon distributions. 



Figure 5b - Two-photon pair production rate versus kT/mc 2 for a thermal 
synchrotron spectrum and different magnetic field strengths. 

Figure 5c - Two-photon pair production rate versus power law index, s, for a 
power law photon spectrum. 

Figure 6a - Ratio of one-photon to two-photon pair production rates at 
temperature kT for a Maxwellian photon distribution. Curves are labeled with 
different values of magnetic field. The dashed line is the thermal 
equilibrium photon density, ngg. 

Figure 6b - Same as Fig. 6a but for a thermal bremsstrahlung spectrum. 

Figure 6c - Same as Fig. 6a but for a thermal synchrotron spectrum. 

Figure 6d - Ratio of one-photon to two-photon pair production rates at 

different magnetic field strengths for a power law spectrum. Curves are 
labeled with power law index. 
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